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Abstract 

We demonstrate with experiments and simulations how microscopic self-propelled particles nav- 
igate through environments presenting complex spatial features, which mimic the conditions inside 
cells, living organisms and future lab-on-a-chip devices. In particular, we show that, in the presence 
of periodic obstacles, microswimmers can steer even perpendicularly to an applied force. Since such 
behaviour is very sensitive to the details of their specific swimming style, it can be employed to 
^ develop advanced sorting, classification and dialysis techniques. 
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Brownian motion is the result of random collisions between a microscopic particle and 
the molecules of the surrounding fluid. In contrast, self-propelled particles additionally take 
up energy from their environment and convert it into directed motion [TH3] ; accordingly, 
their motion is a superposition of random fluctuations and active swimming. Examples of 
such microswimmers range from chemotactic cells |JH6] to artificial systems, where, e.g., 
artificial flagella are put into motion by magnetic fields [71410] or thermophoretic forces [TT] ; 
also, micron-sized Janus particles have been realized where partial coating with a catalyst 
leads to non-isotropic electrochemical reactions and thus to directed motion [T2HH]. Until 
now most studies have concentrated on the behaviour of microswimmers in homogeneous 
environments, where one typically observes a crossover from ballistic motion at short times 
to enhanced diffusion at long times, the latter due to random changes in the swimming 
direction [ToT - TiT] . However, self-propelled particles often move in patterned environments, 
e.g., inside the intestinal tract, which provides the natural habitat of E. coli [5], or during 
bioremediation, where chemotactic bacteria spread through porous polluted soils [TS] . In a 
similar fashion, artificial microswimmers must also reliably perform their tasks in complex 
surroundings, e.g., inside lab-on-a chip devices [19J or living organisms. 

Here, as a first step towards more realistic conditions under which such microswimmers 
are found and will be employed, we study their behaviour in patterned surroundings where 
frequent encounters with obstacles become important. In particular, we investigate their mo- 
tion within environments featuring simple topographical structures such as a straight wall 
and periodically arranged obstacles. Self-propelled Janus particles (Fig. la) are obtained 
from paramagnetic silica spheres with radius R = 2.13 /mi (Microparticles, Si02-MAG- 
S1975). A layer of such particles is deposited onto a substrate. After the solvent evapo- 
rated, their upper sides are coated by thermal evaporation with 2 nm chromium and 20 nm 
gold. This leads to the formation of capped particles. The gold surface is functionalized 
with COOH-terminated thiols (11-Mercapto-undecanoic acid) to render the caps strongly 
hydrophilic. When such Janus particles are suspended in a critical mixture of water and 
2,6-lutidine|20j (0.286 mass per cent of lutidine) below the critical temperature Tc = 307 K, 
they undergo normal Brownian motion (Fig. lb). However, when the entire sample cell is 
homogeneously illuminated with light (A = 532 nm) at low intensities / (< 0.2/iW//im 2 ), 
the particle's motion strongly depends on the incident light intensity as shown in Figs. lc-f. 
In particular the trajectories become ballistic at short times as this is typically found for 
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FIG. 1: (color online) Self-propulsion of Janus particles in critical mixtures, (a) Scanning electron 
microscopy image of a colloidal particle with a 20 nm thick gold cap (highlighted) . (b-f ) 2D trajecto- 
ries (1000s) of a Janus particle for illumination intensities 1 = 0, 69, 92, 115, 161nW//um 2 (respec- 
tively to b, c, d, e and f). (g-i) Probability distribution of the angle between the cap and the particle 
velocity for 1 = (g), 69 (i), and 138nW//xm 2 (i). (j) Experimentally determined mean square 
displacements (symbols) and relative fittings to Eq. (1) (lines) for 1 = 0, 69, 92, 115, 161nW/;um 2 
(respectively squares, circles, up-pointing triangles, down-pointing triangles and diamonds). The 
inset is a magnification of the data for 1 = 0. The solid lines correspond to fits according to Eq. 
(1), see also Tab. 1. 

self-propelled objects. To follow the particles' trajectories with video microscopy, their mo- 
tion is vertically confined by adjusting the sample height to about 7 /mi. Such sample cells 
are made from two microscope cover slides. The porous structures and the walls, which also 
work as spacers, are realized on top of the first slide with soft lithography where a 7/zm 
thick layer of photoresist (MicroChem, SU-8 2005) is first spin-coated and then exposed to 
ultraviolet light through appropriate transmission masks. After development (Micro Resist 
Techn. GmbH, mr Dev 600) and hard backing the structures are filled with colloidal sus- 
pensions and sealed with the second cover slide and epoxy glue. The propulsion is caused by 
diffusiophoresis, which occurs when a particle is subjected to a concentration gradient within 
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FIG. 2: (color online) Encounter between a microswimmer and a wall, (a-e) Time series of snap- 
shots demonstrating the approach (a), contact (b-d) and detachment (e) of a Janus particle and 
a wall. The particle in the left corner is an irreversibly stuck particle which serves as a reference 
position, (f) Experimentally measured (bars) particle-wall contact time distribution p(r c ) between 
a microswimmer and a wall (/ = HSnW/yum 2 ) for more than 250 encounter events. The line 
represents the result of a numerical simulation with L and r taken from the experimental values 
of Tab. 1. 



TABLE I: Characterization of self-propulsion. From the fit of Eq. (1) to the measured MSD we 
obtain the intensity dependent average swimming speed v and time r during which the swimming 
direction is maintained. From v and r we calculate the swimming length L = vt and the effective 
diffusion coefficient D e ff = Dq + j 2 -. The errors correspond to standard deviations. 
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0.031 ± 0.006 
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190 ± 20 
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16 ±3 


0.370 ± 0.069 
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190 ± 20 


175 ± 20 


33 ±8 


1.49 ± 0.42 


115 


220 ± 30 


265 ± 33 


58 ± 14 


3.89 ± 1.01 


138 


240 ± 40 


310 ± 28 


74 ± 13 


5.80 ± 1.49 


161 


230 ± 40 


360 ± 27 


83 ± 12 


7.48 ± 1.86 



the solvent [21]. Here, such gradients are created by the particle itself (self-diffusiophoresis 
[221 [23]) because the incident light is absorbed by the gold cap and thus leads to a local 
demixing of the binary solvent once the cap temperature exceeds Tq- In case of a hydrophilic 
cap this leads to an accumulation of the water-rich phase at the gold cap. 
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From the recorded videos, we also determined the particles time-dependent propulsion di- 
rection v together with their current cap orientation n and the enclosed angle a (see inset 
of Fig. lg). Figs, lg-j show the normalized probability distribution p(a) for increasing 
illumination intensity I. Clearly, p(a) becomes strongly peaked around a = with growing 
/ which demonstrates that the propulsion force acts in the direction opposite to the gold 
cap. When the gold cap is made hydrophobic, the propulsion direction becomes reversed 
since then the lutidine-rich phase accumulates at the cap and thus leads to a change in the 
sign of the concentration gradient. 

The symbols in Fig. lj show the experimentally determined mean square displacements 
(MSD) for the conditions of Figs, lb- f. In the absence of illumination the particles un- 
dergo normal Brownian motion, where the MSD equals Ar 2 = AD^At (inset of Fig.lj) with 
diffusion coefficient Dq = 0.031 ± 0.006 fim 2 /s. With increasing intensity, the MSDs clearly 
deviate from a diffusive behaviour but are well described by 



which generally characterizes the motion of self-propelled particles [151 H7j. Here, r is the 
average timescale over which the trajectory direction is maintained. Accordingly, the aver- 
age velocity of the particle is v = L/t [15] with L the swimming length, i.e. the average 
length of rather straight segments in the particles trajectory. Due to the random changes in 
the particles direction, for long times (At » r) Eq. (1) leads to an effective diffusion with 
D e fj = D + while for short times (At << r) the particles motion becomes ballistic with 
Ar 2 oc L 2 At 2 . The values of the free fitting parameters r and v for which best agreement 
with the data has been achieved are listed together with L and -D e ff Tab.l. From this it 
becomes obvious that L increases (above some threshold value) linear with the illumination 
intensity. In contrast, no significant variation of r with I is observed. The value r « 200 s 
for which best agreement with the experimental data is obtained is close to the timescale 
of rotational diffusion = 1/Dr = 4R 2 /3D Q with Dr the rotational diffusion coefficient. 
This yields tr = 188 s for the R = 2.13/im particles which is indeed close to r. Similar 
agreement (within ±10%) between r and tr has also been found for Janus particle with 
radii between R = 1.0 and 0.5 /im. This suggests that directional changes in the trajectories 
result from cap reorientations due to rotational diffusion. 
When self-propelled particles swim through a patterned medium, in particular for high 
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FIG. 3: (color online) Microswimmers in a patterned environment, (a-e) Typical trajectories 
of self-propelled particles moving through a triangular lattice (lattice constant L c = 35 /tm) of 
elliptical obstacles when a drift force F = 0.12 pN is applied along the y-direction. (a) Brownian 
particle (no propulsion), (b) L = 16 /im, (c) L = 24 /mi, (d) L = 33 /xm and (e) L = 83 /mi. 
(f-j) Corresponding histograms of the experimentally measured (red/light grey) and simulated 
(blue/dark grey) directions of the particle trajectories as defined by two points in the trajectory 
separated by 100 /mi. The parameters for the simulations are taken from Tab. 1. 

swimming velocities, frequent encounters with obstacles will occur. Figs. 2a-e show the 
interaction between a Janus particle and a wall. First the particle approaches the wall (Fig. 
2a) and gets in contact (Fig. 2b). Then it slides along the wall (Fig. 2c) until rotational 
diffusion realigns the particle so that its orientation vector n points away from the wall and 
leads to particle detachment (Figs. 2d,e). Measuring the distribution of the particle- wall 
contact time r c , we find a monotonic decrease as shown in Fig. 2f (bars). We also performed 
numerical Brownian dynamics simulations [23], where the motion of the Janus particle is 
modelled by a superposition of random diffusion and ballistic motion with velocity v as 
inferred from Tab.l. The direction of the ballistic component is determined by the orienta- 
tion of the particle, which is continuously realigned due to rotational diffusion. Whenever 
the particle encounters a wall, the translational motion perpendicular to the wall is set to 
zero, while the horizontal component remains unchanged. From the simulated particle tra- 
jectories we calculated the corresponding distribution p(r c ) (line in Fig. 2f), which shows 
good agreement with the experimental data. This suggests that the particle-wall encounter 
mechanism is correctly described by this simple model and, in particular, that the rotational 
diffusion remains largely unaffected by the proximity to the wall. 

We have also investigated the behaviour of self-propelled particles in the presence of a two- 
dimensional periodic pattern where straight unlimited swims are only possible along certain 
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directions. We have chosen a structure made of a series of ellipsoidal pillars arranged in 
a triangular lattice (lattice constant L c = 35 /im, Fig. 3a). Within such structures, long 
swimming cycles are only possible along two main directions: at ±60° and ±90° with respect 
to the y-axis. Otherwise the motion is strongly hindered due to collisions with the obstacles. 
In the presence of an additionally applied drift force, this leads to strong differences in the 
particle trajectories depending on their swimming length. In our experiments a constant 
drift force acting on the paramagnetic particles in the y-direction is generated by a magnetic 
field gradient. From the average drift speed of a Brownian particle (/ = 0/iW//xm 2 ), the 
drift velocity has been determined to Vd = 0.97 /mi/s, which corresponds to a Stokes force 
of F = 0.12 pN. 

The typical trajectory of a Brownian particle is shown in Fig. 3a. Because of the Peclet 
number Pe « 1000, the effect of the diffusion is rather weak and the particle meanders al- 
most deterministically through the structure in the direction of F. For increasing swimming 
lengths, however, significant changes in the shape of the trajectories are observed. This 
becomes particularly pronounced for L > L c , where the particles perform swimming cycles 
of increasing length along the diagonal channels (Figs. 3c, d). For L = 83 /im the propulsion 
becomes so strong that the particles partially move perpendicular to the drift force (Fig. 
3e); occasionally even motion against the drift force can be observed. 

The direction of the particle motion through the structure is characterized by the direction 
(with respect to the y-axis) of the line connecting points of the trajectory separated by a 
distance of 100 /im. The probability distributions of these angles are shown by the red/light 
grey polar histograms in Figs. 3f-j. One clearly observes that with increasing L the propaga- 
tion of particles along the direction of the applied drift becomes less likely, while trajectories 
along ±60°, i.e. along the directions that permit long swimming events, become more fre- 
quent. We remark that, differently from the deflection of Brownian particles in a periodic 
potential [2l)H2"7] . this mechanism relies on the dynamical properties of the micros wimmers. 
We also compare these results with numerical simulations (blue/dark grey polar histograms 
in Figs. 3f-j), which show good agreement with the experimental data. 

With the additional possibility of varying the drift force, these observations can be ex- 
ploited to spatially separate self-propelled particles with small differences in their individual 
swimming behaviour. This is demonstrated in Fig. 4, where we show the deflection effi- 
ciency as a function of L, as defined by the probability that the mean particle trajectory is 
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FIG. 4: (color online) Deflection efficiency. Measured probability that particles are deflected by 
more than 30° after a travelling length of 100 fim as a function of the swimming length L for 
various imposed magnetic drift forces F = 0.06 db 0.02 pN (diamonds), 0.12 ± 0.05 pN (squares) and 
0.28±0.12pN (triangles). The solid lines are the results of numerical calculations. The parameters 
for the simulations are taken from Tab. 1. 

deflected by more than 30° after a travelling length of 100 fim. As symbols (lines) we have 
plotted experimental (simulation) data obtained for different drift forces. The black data 
correspond to F = 0.12 pN (Fig. 3). Here the deflection efficiency shows a strong increase 
around L w 30/zm and a flattening towards larger swimming lengths. The blue and red 
data are obtained for F = 0.06 and 0.28 pN, respectively, and demonstrate that the sorting 
efficiency strongly responds to variations in F. Accordingly, the deflection (transmission) 
of self-propelled particles while crossing a patterned structure can be easily tuned by the 
appropriate choice of F. 

Finally, we remark that the sorting mechanism discussed here can be directly applied to 
other self-propelled objects. In these cases, drift forces can be created, e.g., by electric fields 
or by a solvent flow through the device. Therefore, we expect that our method may find 
wider use as an efficient technique to characterize the behaviour of motile cells and bacteria 
even when the details of their swimming behaviour differs from those of the particles em- 
ployes in this work. Compared to other concepts for the sorting of chemotactic bacteria, our 
strategy avoids the need for the creation of chemical gradients which are difficult to stabilize 
in time and space ES] ■ 
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